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Resumen 


Sabat, 1, Link, O., £ Ettmer, B. (enero-febrero, 2015). Flujo 
inducido por el rompimiento de una presa: análisis mediante 
modelación numérica. Tecnología y Ciencias del Agua, 6(1), 
5-23. 


Se aplica el método de los volúmenes finitos para dar 
solución a las ecuaciones de Saint Venant 2D con términos 
fuente que gobiernan el escurrimiento de aguas someras, a 
fin de calcular el flujo inducido por el rompimiento de una 
presa. El modelo se verifica mediante el cálculo de siete 
casos con solución conocida, que son: (1) establecimiento 
del flujo normal a partir de una superficie libre horizontal; 
(2) establecimiento del flujo gradualmente variado en un 
canal con cambio de rugosidad; (3) escurrimiento sobre un 
obstáculo; (4) formación del salto hidráulico oblicuo; (5) 
oscilación amortiguada de la superficie libre en un estanque 
parabólico; (6) escurrimiento inducido por el rompimiento 
de una presa circular, y (7) escurrimiento inducido por el 
rompimiento parcial de una presa. El modelo se aplica para 
simular el flujo inducido por el rompimiento instantáneo 
de presas, a fin de analizar el efecto de la configuración 
geométrica del conjunto presa, embalse, valle y cauce 
sobre el hidrograma de salida. Se simularon 28 escenarios, 
considerando un valle con sección parabólica y un cauce con 
sección trapezoidal, donde se varió la altura y orientación 
del muro de la presa, el ancho de la brecha, la pendiente 
longitudinal y la rugosidad del cauce aguas abajo. Los 
resultados obtenidos muestran que la altura de la presa es la 
variable más influyente en el caudal máximo del hidrograma 
de salida y que el ancho de la brecha tiene un efecto 
importante en rompimientos parciales. La rugosidad del 
lecho atenúa de forma considerable el caudal máximo hacia 
aguas abajo y aumenta el tiempo de crecida. La orientación 
del muro y la pendiente longitudinal tienen un efecto 
menor sobre las crecidas provocadas por rompimiento. 
El modelo desarrollado predice de manera adecuada el 
complejo comportamiento que exhibe el flujo inducido por 
el rompimiento instantáneo de una presa, constituyendo una 


herramienta de simulación con capacidad predictiva. 


Palabras clave: rompimiento de presas, ecuaciones de Saint 
Venant 2D, hidráulica computacional, esquemas de shock. 


Abstract 


Sabat, I., Link, O., € Ettmer, B. (January-February, 2015). Flow 
Caused by Dam Breaks: Numerical Modeling Analysis. Water 
Technology and Sciences (in Spanish), 6(1), 5-23. 


The finite volumes method is applied to solve 2D Saint Venant 
equations with source terms that determine runoff of shallow 
water, in order to calculate flows caused by dam breaks. The 
model is verified by calculating six cases with known solutions: 
(1) establishing normal flow based on a horizontal free surface; (2) 
establishing flow that is gradually varied in a channel with changes 
in roughness; (3) runoff over an obstacle; (4) formation ofan oblique 
hydraulic jump; (5) damped oscillation of a free surface in a parabolic 
section; (6) runoff from the breaking of a circular dam and (7) runoff 
from the partial breaking of a dam. The model is applied to simulate 
the flow caused by the instantaneous breaking of dams in order to 
analyze the effect of the joint geometric configuration of the dam, 
reservoir, valley and channel on the outflow hydrograph. A total 
of 28 scenarios were simulated for a valley with a parabolic section 
and a channel with a trapezoidal section, in which the height and 
position ofthe dam wall were varied as well as the gap width, length 
of the slope and roughness of the channel downstream. The results 
obtained show that the height of the dam is the variable that has the 
largest effect on the peak flow of the outflow hydrograph and the gap 
width has a large effect on partial breaks. The roughness of the bed 
considerably mitigates the peak flow downstream and increases flood 
time. The position of the mural and the length of the slope have a 
smaller effect on floods caused by dam breaks. The model developed 
adequately predicts the complex behavior of the flow caused by the 
instant breaking of a dam and is therefore a simulation tool with 
predictive capabilities. 


Keywords: Dam break, Saint Venant 2D equations, computational 
hydraulics, shock capturing conditions. 
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Introducción 


El rompimiento de presas es un tema relevante 
en ingeniería hidráulica debido a la alta 
probabilidad de ocurrencia de consecuencias 
catastróficas que provoca el flujo inducido 
aguas abajo del muro. Singh (1996) registró 
cerca de mil casos de rompimientos de presa 
sucedidos desde el siglo XII, de los cuales cerca 
de 200 se registraron en el siglo XX. Muchos 
rompimientos se encuentran bien docu- 
mentados en la literatura especializada, como 
el de las presas Malpasset (Hervouet é Petitjean, 
1999), St. Francis (Begnudelli $: Sanders, 2007), 
Tous (Alcrudo $ Mulet, 2007) y Gleno (Pilotti, 
Maranzoni, Tomirotti, £ Valerio, 2011). Singh 
(1996), y Vischer y Hager (1998) indican que las 
principales causas de rompimiento de presas 
han sido problemas en las cimentaciones (ca., 
37%) y overtopping (ca., 30%). Dada la relevancia 
del problema, muchos países exigen hoy día 
como parte del proyecto de un embalse, el 
análisis de seguridad de la presa, que incluye 
el cálculo del flujo inducido por su eventual 
rompimiento. 

De manera típica, las presas de tierra y 
enrocado presentan rompimientos graduales, 
con una abertura de brecha que demora 
del orden de horas, y que induce un flujo 
gradualmente variado hacia aguas abajo, 
mientras que las presas de hormigón presentan 
rompimientos instantáneos, con la formación 
de una onda negativa hacia aguas arriba, 
controlada por la geometría del valle, y una 
onda de shock que se propaga hacia aguas 
abajo (e.g., Vischer e Hager, 1998). 

Rubio, Martínez y Meza (2012) estimaron 
el caudal máximo causado por el rompimiento 
gradual de una presa aplicando distintas 
formulaciones empíricas. Realizaron un rastreo 
de la crecida hacia aguas abajo utilizando el 
software HEC-RAS, considerando un canal 
trapezoidal, a fin de establecer la distancia a 
la que el paso de la avenida ya no representa 
peligro para la población. 

La condición más desfavorable es el rom- 
pimiento instantáneo, que puede modelarse 
mediante las ecuaciones de Saint-Venant, inclu- 


yendo los términos fuente que representan el 
efecto de la pendiente del cauce y de la fricción 
en el lecho sobre el escurrimiento. Las solucio- 
nes analíticas conocidas se limitan al caso 1D 
para un embalse semi-infinito y canal rectan- 
gular aguas abajo de la presa. Dressler (1952, 
1954) y Whitham (1955) propusieron solucio- 
nes para el rompimiento instantáneo y un ca- 
nal seco con pendiente horizontal. Hunt (1982, 
1984) resolvió el rompimiento instantáneo para 
un canal seco con pendiente longitudinal; de 
forma reciente, Chanson (2009) propuso solu- 
ciones aplicando el método de las característi- 
cas, para el rompimiento semi-instantáneo de 
un embalse semi-infinito, con canal rectangular 
inicialmente seco. 

En casos reales, la geometría y resistencia al 
escurrimiento del valle y del cauce controlan 
la propagación de las ondas hacia aguas arriba 
y hacia aguas abajo. Por ello, en la práctica 
es pertinente una representación 2DH del 
fenómeno, que se realiza mediante aplicación 
de métodos numéricos. La existencia de ondas 
de shock requiere implementar esquemas de 
discretización espacial con muy baja difusión 
numérica, que garanticen conservación de los 
flujos a nivel local y estabilidad de la solución, 
en especial en la vecindad del shock. La 
mayoría de estos esquemas se ha propuesto 
para el problema de Riemann que no considera 
los términos fuente (e.2., Fraccarollo € Toro, 
1995; Brufau € García-Navarro, 2003). 

Para la solución de las ecuaciones de Saint- 
Venant con términos fuente en presencia de 
ondas de shock, en los últimos años se han 
desarrollado esquemas de discretización es- 
pecíficos, llamados esquemas de shock, los 
que, combinados con una adecuada estrategia 
secado-mojado, presentan una excelente capa- 
cidad de simulación del flujo inducido por el 
rompimiento de una presa. Gottardi y Venu- 
telli (2004) extendieron la formulación del es- 
quema centrado de Kurganov y Tadmor (2000), 
para la discretización de los términos fuente 
en las ecuaciones de Saint Venant, presentan- 
do resultados favorables en la simulación del 
rompimiento de presas idealizadas. Zhu, Vis- 
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ser y Vrijling (2004) implementaron un méto- 
do predictor-corrector basado en el esquema 
MUSCLE-Hancock (van Leer, 1985) para la 
simulación de casos de estudio desarrollados 
en la iniciativa CADAM (Morris, 2001). Soares- 
Frazáo y Zech (2008) presentaron simulaciones 
numéricas de los efectos del flujo inducido por 
un rompimiento sobre una ciudad idealizada 
compuesta por una grilla de 5 x 5 edificios cua- 
drados, que se orientó a 45 y 90” respecto a la 
dirección principal de la corriente. Los cálculos 
se compararon con resultados experimentales, 
obteniendo una buena concordancia entre en- 
sayos de laboratorio y resultados numéricos. 
Wang, Liang, Kesserwani y Hall (2011) exten- 
dieron la formulación 1D propuesta por Liang 
y Marche (2009), para simular rompimientos 
idealizados y también el de la presa Malpasset. 

En este artículo se extiende el modelo nu- 
mérico desarrollado por Link y Donoso (2008) 
al caso 2D con términos fuente y estrategia 
secado-mojado. El modelo se verifica mediante 
el cálculo de siete casos con solución conocida, 
que son: 1) establecimiento del flujo normal, a 
partir de una superficie libre horizontal; 2) es- 
tablecimiento del flujo gradualmente variado 
en un canal con cambio de rugosidad; 3) escu- 
rrimiento sobre un obstáculo; 4) formación del 
salto hidráulico oblicuo; 5) oscilación amor- 
tiguada de la superficie libre en un estanque 
parabólico; 6) escurrimiento inducido por el 
rompimiento de una presa circular, y (7) escu- 
rrimiento inducido por el rompimiento parcial 
de una presa. El modelo se aplica para simular 
el flujo inducido por el rompimiento instantá- 
neo de una presa, a fin de analizar el efecto de 
la geometría del conjunto presa, embalse, valle 
y cauce sobre el hidrograma de salida. 


Metodología 
Ecuaciones gobernantes 
El sistema de ecuaciones gobernantes corres- 


ponde a las ecuaciones de Saint-Venant con 
términos fuente, que en forma vectorial es: 


W iv.ps (1) 
ot 
aw, oE(W) ñ 0G(W)_ A 0) 
ot Ox dy 


donde W es el vector de variables de flujo; F y G 
son vectores de flujo que componen la matriz de 
flujo F en dirección x e y, respectivamente, y S es 
el vector de términos fuente correspondientes a 
la pendiente del lecho y la fricción del fondo 
del cauce. Los vectores W, F, G y S están dados 
por: 


, u:h , 
-h 
W=| uh | F=| 4" -h+£ ' 
v:h 2 
u:v:h 
v:h 
G= u:v:h ] 
.y? 
+8 h 
2 
0 
S=  g:h(S,,-S;) (3) 
8 (Su fy 


donde h (m) es la profundidad de la columna 
de agua; u y v (m/s), la velocidad integrada 
en la vertical en las direcciones de los ejes 
cartesianos x, y; t (s), el tiempo; g (ms?) 
aceleración de gravedad; x (m), la distancia en 
sentido longitudinal; y (m) es la distancia en 
sentido transversal. Las pendientes de fondo 
Y y (m/m) son: 


(4) 


donde z, corresponde a la cota batimétrica. La 
pendiente friccional S,, y S,, (m/m) se calcula a 
partir de la ecuación de Manning como: 
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4 
hH3 
ven du? + v 
2 ———— A (5) 
3 


donde 1 (m1? s) es el coeficiente de rugosidad 
de Mamning. 


Formulación discreta del problema 


Las ecuaciones gobernantes se llevaron a la 
forma discreta aplicando el método de los 
volúmenes finitos, utilizando un esquema 
descentrado perteneciente a la familia de 
Q-esquemas, introducida por Harten, Lax y Van 
Leer (1983), y extendida para las ecuaciones 
de aguas someras con términos fuente por 
Bermúdez y Vázquez (1994) y Bermúdez, 
Dervieux,  Desiderii (1998). 
Integrando la ecuación gobernante se tiene: 


y Vázquez 


aW 
Jr joas (vr (waa - $(8)os (6) 
a(27)eJ05=5:a (7) 
dt) % 
IWY. luz 
d )--¿3hises (8) 


donde A es el área; S, el contorno del volumen 
de control, y n es un vector normal al contorno. 
El término de flujo se aproximó mediante la 
expresión: 


JEnas- Sr, -o(W! Win) 0) 


j=l 


donde T ; (m) es la longitud de la cara entre los 
volúmenes i y j, y la función Y se define como: 


p(U,V,n) z a 
1 
=> ¡e(Wa(U,V),n)|.(V-u) (10) 
donde: 

Z(U,n)=A'(U,n)-U (11) 
wa(uv)- (2%) (12) 
Q(W,n) =A (W,n) (13) 
A (W,n)=X-A-X” (14) 


U y V son vectores de propiedades del flujo de 
volúmenes finitos; n es el vector normal a la 
cara entre U y V, y A* es la matriz Jacobiana, 
según Vijayasundaram (1986): 


0 1 1 
X=| n, u+cn, u-c"n, (15) 
=N, U+CN, UCM, 
A .0 0 
A=| 0 A 0 (16) 
0.04 


Nótese que A,, A y A,(ms*) corresponden 
a los autovalores de A*. Para evitar la 
indeterminación del problema en las discon- 
tinuidades, se introdujo la corrección de 


entropía de Harten (1984): 


0.5-sign(a,): Bl. silk;|se 
. le 3-0)07) 


» sila.| >e 


1 
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Se consideró e = 10% (m s”). 
Términos fuente 


La pendiente de fondo se discretizó como: 


a "EA, Z ¡72 ¡ 
Sor == a > 2 A, iS d, ] cos(8,), 
0 MERA Zo o. 
Soy = > a 2 A, a d, ] ssin(9,) 50) 


donde z,, y 2, Corresponden a la cota batimétrica 
del volumen ¡ y su vecino j, respectivamente; 
d, es la distancia entre los centros de los 
volúmenes; 0,, el ángulo generado por la 
horizontal y la línea que une los centros de los 
volúmenes i y j, medido en sentido antihorario; 
A,, el área del triángulo formado por el centro 
del volumen ¡ y los vértices comunes entre 
los volúmenes ¡ y j, y A, es la superficie del 
volumen i. La pendiente friccional se discretizó 
mediante la ecuación de Manning, como: 


Un; ca fu; +0; 
Spj= 4 / 
h3 
2 2 2 
DN; “qu, +0; 
Sy NA (19) 


donde u, y v,sonlas componentes de la velocidad 
en la dirección del eje x e y, respectivamente; 
h, es la profundidad de la columna de agua, y 
n, corresponde al coeficiente de rugosidad de 
Mamning en el centro del volumen . 

Las derivadas temporales se aproximaron 
utilizando el esquema explícito de Euler de 
primer orden, de forma que la solución discreta 
es: 


r,0(91,01) 


0 
Meecinos AZ ZO 
j 0 0 
+9*h, —— L-cos(8,,) 
j A, d, 
Moecinos A. o 
AAN -sin(8,) 
A J 
] 1 y 
0 
2 
UN; ca fu; +0 
a 4 
h3 
1 
0 Mica fui RO 
4 
h3 (20) 


Condiciones de borde 


Para el borde abierto se impuso el contorno 
transmisivo mediante un vecino ficticio con 
propiedades del flujo idénticas a las del 
volumen ubicado en el contorno (Toro, 2001): 


W",=W, (21) 


donde W”, corresponde al elemento ficticio 
y W, al volumen en la frontera del dominio. 
Para el borde cerrado se impuso el contorno 
reflectivo, donde la velocidad en la dirección 
perpendicular a la cara se invierte (Toro, 2001): 


ll cos(0,,)  sin(0,) | a, 


-sin(8,,) cos(0,.) 0, | (23) 
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(24) 


donde h, u, y v, corresponden a las propiedades 
del volumen en la frontera del dominio; 0, es el 
ángulo perpendicular a la cara en la frontera y 
saliente del volumen; u”, y v”, son velocidades 
rotadas en 0,, y u”, y v”,son las componentes 
de la velocidad del vecino ficticio en los ejes 
cartesianos. 


Estrategia secado/mojado 


La estrategia secado /mojadoimponela nulidad 
de los Jacobianos asociados con un volumen 
seco cuando existen dos vecinos: uno seco y 
uno mojado. En este caso, las profundidades 
mínimas se limitan a valores mayores o iguales 
a cero, y se realiza un balance de masa entre 
los vecinos mojados, a fin de garantizar la 
continuidad. Matemáticamente: 


0, Vh,<0 
M,= Ñ (25) 
| L Vh,>0 
3, = min(0,h,) (26) 
M,:A,*8, 
hy =h + (27) 
My i 
] 
h =h, - 0, (28) 
h 
W, = i 1 (29) 
v. 


donde M, indica si el volumen ¡ está seco o 
mojado y d,es la cantidad de masa a balancear. 
La figura 1 muestra el diagrama de flujo del 
modelo propuesto. 


Resultados 
Verificación del modelo 


La capacidad del modelo desarrollado se 
verificó mediante siete benchmark test: 1) 
establecimiento del flujo normal a partir de una 
superficie libre horizontal; 2) establecimiento 
del flujo gradualmente variado en un canal con 
cambio de rugosidad; 3) escurrimiento sobre 
un obstáculo; 4) formación del salto hidráulico 
oblicuo; 5) oscilación amortiguada de la 
superficie libre en un estanque parabólico; 6) 
escurrimiento inducido por el rompimiento de 
una presa circular, y 7) escurrimiento inducido 
por el rompimiento parcial de una presa. 
Los siete casos mencionados tienen solución 
conocida y permiten evaluar distintos aspectos 
del modelo numérico, como la convergencia a 
estados de equilibrio, conservación, propiedad 
C de Bermúdez y Vázquez (1994), simulación 
de discontinuidades en el flujo y modelación 
del proceso de secado-mojado. Los resultados 
numéricos se compararon con soluciones 
analíticas, semi-analíticas o empíricas en forma 
gráfica, e indicando la diferencia porcentual 
máxima. 


Establecimiento del flujo normal a partir de 
una superficie libre horizontal 


Se realizó la simulación del escurrimiento en un 
canal rectangular de ancho unitario constante 
de largo igual a 100 m, pendiente igual a 1/2 000 
m/m y un coeficiente de Manning igual a 
0.02 sm”. Las condiciones inicial y de borde 
fueron: 

h 


=0.5-0.02:(100-x) q =0 


(x,y,0) x(x,y,0) AS Uylxy.0) 


Moya 205 yg y) 2005" (ao)? Uxli-ax y ¡a 


La solución está dada por la ecuación de 
Mamning y es igual a 0.22 m*/s/m. Se utilizó 
una discretización espacial Ax = 0.25 m y Ay 
= 0.25 m, y un paso temporal At = 0.05 s, con 
un número de Courant CFL = 0.15. La figura 
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Ingreso de datos: dominio, cálculo, malla, 
rugosidad, tiempo de simulación, At 


Cálculo de la geometría de la malla 
Impone condición inicial 
Inicio ciclo temporal 0: At: 


Inicio ciclo espacial 


Cálculo del flujo en las caras de cada volumen 
Cálculo de términos fuentes en cada volumen 


Impone condiciones de borde 
Estrategia de secado / mojado 


in ciclo temporal 
Guarda resultados 


Fin ciclo espacial 
Gráfica variables u, v y h 


Tm > 


Figura 1. Diagrama de flujo. 


2 muestra la malla, condición inicial y eje 
hidráulico calculado. 

El cuadro 1 muestra la diferencia porcentual 
máxima entre la profundidad y caudales 
calculados por el modelo y la ecuación de 
Mamning. Se observa que la predicción del 
modelo prácticamente coincide con la de 
Mamning. 


Establecimiento del flujo gradualmente variado en 
un canal con cambio de rugosidad 


Se realizó la simulación del escurrimiento 
gradualmente variado que causa un cambio 
de rugosidad en un canal rectangular de 
ancho unitario constante, de largo igual a 
100 m y pendiente igual a 2/1000 m/m. El 
canal se dividió en tres tramos de 33.3 m, con 
rugosidades de Manning iguales a 0.02, 0.04 y 
0.02 sm", respectivamente. Las condiciones 
inicial y de borde fueron: 


h,..., y =0.92-0.002-(100- x) Oca) Ia) O 


(4.40) 


Ixoye 2921 Gary =05* (a «(Ly zas) Y ista] 


h =0.6 


(L.y) 


La solución numérica se comparó con 
aquella calculada por el software HEC-RAS, 
que resuelve las ecuaciones de Bresse mediante 
el método del paso estándar. Se utilizó una 
discretización espacial Ax = 1 m y Ay = 0.4 m, 
y un paso temporal Af = 0.10 s, con un número 
de Courant CFL = 0.14. La figura 3 muestra 
la malla, condición inicial y eje hidráulico 
calculado. 

Los resultados de la simulación muestran un 
eje hidráulico semejante al obtenido mediante 
HEC-RAS. El cuadro 1 muestra la diferencia 
porcentual máxima entre la profundidad y 
caudales calculados por el modelo y el HEC- 
RAS, que no alcanza el 0.5%. 


Escurrimiento sobre un obstáculo 


Se realizó la simulación del escurrimiento sobre 
un obstáculo con forma de parábola, ubicado 
en el fondo de un canal rectangular con ancho 
unitario de 25 m largo, sin pendiente de fondo 
ni rugosidad. La geometría del obstáculo fue: 


O mparax<7.5m 
2 
-10 
z=30.2 a) m para7.5m<x<12.5 m Vy 


Om para x > 12.5 m 
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0 


0 10 20 30 40 50 60 70 80 90 100 


x (m) 


== Ecuación de Manning 


Modelo propuesto 


10 20 30 40 50 60 70 80 90 100 
x (m) 


0 


Figura 2. Malla y condición inicial (izquierda) y eje hidráulico calculado (derecha) con el modelo y la ecuación de Manning. 


Cuadro 1. Error de los resultados del modelo con respecto a las soluciones conocidas en el cálculo de la profundidad y caudal. 


de Error promedio (%) Error máximo (%) 
h q h q 
Test 1 0.01 0.00 0.03 0.00 
Test 2 0.36 0.02 0.76 0.04 
Test 3a 0.05 0.04 0.33 0.19 
Test 3b 1.98 0.03 10.07 0.27 
Test 3c 2.54 0.03 8.80 0.35 
Test 4 0.01 N/A 1.12 N/A 
Test 5 1.18 1.37* 3.20 642.16* 


* Error asociado con la velocidad del flujo. 
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Figura 3. Malla y condición inicial (izquierda) y eje hidráulico calculado (derecha) con el modelo y con el software HEC-RAS. 


Se impusieron tres condiciones iniciales y 
de borde distintas, a fin de simular: a) el flujo 
subcrítico sobre un obstáculo; b) el cambio de 
régimen sub a supercrítico, y c) la formación 
del salto hidráulico, mostradas en el cuadro 2. 

Los resultados numéricos se compararon 
con la solución analítica publicada por 
Alcrudo y Benkhaldoun (2001). Se utilizó una 


discretización espacial Ax = 0.1 m y Ay=0.1 m, 
y un paso temporal At =0. 01 s, con un número 
de Courant CFL = 0.26, 0.32 y 0.38 en los casos 
a, b y c, respectivamente. La figura 4 muestra 
la condición inicial y los ejes hidráulicos 
calculados. 

Los resultados de la simulación muestran 
un eje hidráulico semejante al obtenido me- 
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Cuadro 2. Condiciones iniciales y de borde para los test 3a, 3b y 3c. 


Test Condiciones iniciales (CI) Condiciones de borde (CB) 
Mio) =2.5 UN (0,91) = 4.42, Uy(0,1,1) =0 
Test 3a 
UR x,y,0) = Uy(x,10) =0 Mo) =20 
M0) =1.25 eo) =1.53, yo) =0 
Test 3b 
UR 140) = Uy(x:10) =0 vo = Was a + | 
Ko) =0.45 Ix(o0y:) 0.18, Iy(0,y:) 0 
Test 3c 
Axa) = Tyco) 0 oy = 0.33 


diante Alcrudo y Benkhaldoun (2001). El cua- 
dro 1 muestra la diferencia porcentual máxima 
entre la profundidad y caudales calculados por 
el modelo, y la solución de Alcrudo y Benkhal- 
doun (2001). En el caso a), los errores son muy 
pequeños, inferiores al 0.5%. En los casos b) y 
c), la profundidad calculada presenta errores 
de hasta un 10.87%, lo que se atribuye a las 
presiones dinámicas que se generan sobre el 
obstáculo y que el modelo asume como hidros- 
táticas. 


Formación del salto hidráulico oblicuo 


Se realizó la simulación del escurrimiento en 
un canal rectangular con 60 m de longitud, 40 
m de ancho, sin pendiente ni roce, con una con- 
tracción que se inicia 20 m aguas abajo de la en- 
trada, con un ángulo de 8.95”. Las condiciones 
iniciales y de borde fueron: 


Condiciones iniciales 


h =1 


(+40) 


UN (x,4,0) = 8.57 9, (x,y,0) a 0 


Condiciones de borde 


h =1 


(0,y,1) 


Ix(0,y,) = 8.57 Uy (0,4) 5 0 


y garantizan un escurrimiento supercrítico 
antes del angostamiento. La solución numérica 
se comparó con la solución analítica publicada 
por Alcrudo y García-Navarro (1993). Se utilizó 
una discretización espacial Ax = Ay = 0.5 m y 
un paso temporal At = 0.025 s, con un número 
de Courant CFL = 0.43. En total, se simularon 
150 s. La figura 5 muestra la malla de cálculo, 
condición inicial y superficie libre calculada. 


SS 
10 15 
x (m) 


x (m) 


x (m) 


Figura 4. Condición inicial (arriba) y ejes hidráulicos calculados para test 3a (izquierda, abajo); test 3b (centro, abajo), 


y test 3c (derecha, abajo), con el modelo y la solución analítica de Alcrudo y Benkhaldoun (2001). 
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La solución analítica entrega alturas conju- 
gadas iguales a 1.0 y 1.5 m, formando un 
ángulo de 30” con respecto al eje longitudinal. 
Se observa que los resultados obtenidos con 
el modelo propuesto coinciden con dicha 
solución. El cuadro 1 muestra la diferencia 
porcentual máxima entre la profundidad 
calculada por el modelo y la solución analítica 
de Alcrudo y García-Navarro (1993). El caudal 
no se compara, ya que éste no se conoce a partir 
de la solución analítica. Puede notarse que las 
diferencias no superan el 1.2%. 


Oscilación de la superficie libre en un recipiente 
parabólico 


Se simuló la oscilación de la superficie libre en 
un recipiente parabólico de radio superficial 
igual a 4 500 m. Inicialmente, se impuso una 
superficie libre inclinada, según Wang et al. 
(2011): 


=h l B?p" a 1 pg, u/2 


Poy) 0 28 g 
(sins) + scos(s)] -Xx- pl 


ES 
(eos(s) E ssin(s) y 


Be" 


Uy = Be" «sin/(st) Uy) = pere -cos(st) 
donde: 
2: 2 
5 pot de 82 o 
2 a 


Se utilizó h, = 10 m, 1=3000 m, B=5 m/s 
y T= 0.002. La solución numérica se comparó 
con la solución analítica publicada por Wang et 
al. (2011). Se utilizó una discretización espacial 
Ax = Ay = 50 m, y un paso temporal At = 0.25 s, 
con un número de Courant CFL= 0.03. En total, 
se simularon 3 600 s, que abarca tres periodos 
de la oscilación. La figura 6 muestra la malla 
de cálculo y condición inicial, velograma en el 
punto (x, y) =(1 000,0) para ambas componentes 
u y v, así como la superficie libre calculada para 
t= 900, 1 800, 2 700 y 3 600 s. 

Se observa que el modelo simula de manera 
adecuada las oscilaciones amortiguadas por 
efecto de la fricción, y que los resultados 
coinciden muy bien con la solución analítica de 
Wang et al. (2011). El velograma calculado para 
la componente longitudinal de la velocidad 
coincide con la solución analítica, mientras que 
para la componente transversal presenta un 
desfase con la solución analítica que eleva el 
error máximo. El cuadro 1 muestra la diferencia 
porcentual máxima entre la profundidad 
y velocidad calculados por el modelo, y la 
solución analítica de Wang et al. (2011). 


h(m) 


x (m) 


h (m) 


Figura 5. Malla y condición inicial (izquierda) y superficie libre calculada (derecha) con el modelo. 
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Figura 6. Malla y condición inicial (arriba, izquierda), velograma en el punto (x, y) = (1 000.0) (arriba, derecha) 
y superficie libre calculada para t = 900, 1 800, 2 700 y 3 600 s calculada con el modelo propuesto y la solución analítica 
de Wang et al. (2011). 


Escurrimiento inducido por el rompimiento de una 
presa circular 


Se realizó la simulación del escurrimiento 
inducido por el rompimiento instantáneo de 
una presa con 10 m de altura y forma anular. 
Alrededor de la presa se consideró una 
planicie de inundación horizontal sin fricción, 
con profundidad inicial igual a 1.0 m. Las 
condiciones iniciales y de borde se muestran 
en el cuadro 3. 


La solución numérica se comparó en 
forma gráfica con los resultados de Mingham 
y Causon (1998). No existen antecedentes 
suficientes para comparar valores puntuales 
ni hacer el cálculo de errores asociados con 
las simulaciones. Se utilizó una discretización 
espacial Ax = Ay = 0.25 m y un paso temporal 
At=0.01 s. En total, se simularon 300 s. La figura 
7 muestra la malla de cálculo, condición inicial 
y superficie libre calculada. Los resultados de la 
simulación coinciden con la solución numérica 
de Mingham y Causon (1998). 
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Cuadro 3. Condiciones iniciales y de borde. 


Condiciones iniciales 


10, yx*+y?<=11 


1 y“i+y >11 


Ue(x,y0) = Vyfxy oy 0 


(13.0) — 


Condiciones de borde 


Ax (Borde, 1) = A x(Borde, 1) 


2 y (Borde, e) “Ay (Borde...) 
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Figura 7. Malla y condición inicial (arriba, izquierda), campo de velocidades (arriba, derecha) y superficie libre para t = 0.69 


(abajo, izquierda) y 300 s (abajo, derecha), calculada con el modelo propuesto. 


Escurrimiento inducido por el rompimiento parcial 
de una presa 


Se simuló el escurrimiento inducido por el 
rompimiento instantáneo y asimétrico de una 
presa orientada en ángulo perpendicular al eje 
principal del valle, con 200 m de largo, 10 m 
de ancho y 11 m de altura. La zona de rotura 
se ubicó entre los 95 y 170 m. Aguas abajo de 


la presa se modeló una planicie de inundación 
horizontal sin fricción en dos situaciones: a) 
inicialmente seca, y b) inicialmente con una 
profundidad de agua de 5 m. Las condiciones 
inicial y de borde se muestran en el cuadro 4. 
La solución numérica se comparó en 
forma gráfica con los resultados de Mingham 
y Causon (1998). No existen antecedentes 
suficientes para realizar una comparación 
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de valores puntuales ni para llevar a cabo el 
cálculo del error asociado con las simulaciones. 
Se utilizó una discretización espacial Ax = 5 
m y Ay = 5 m, y un paso temporal Af = 0.01 
s. En total se simularon 150 s. La figura 8 
muestra la malla, condición inicial, superficie 
libre calculada, isolíneas de profundidad y 
campo de velocidades para las situaciones 
a) y b), respectivamente. Los resultados de la 
simulación coinciden con la solución numérica 
de Mingham y Causon (1998). 


Análisis del flujo inducido por el 
rompimiento de una presa 


Se analizó el flujo inducido por el rompimiento 
instantáneo de una presa en un sistema 
hidráulico idealizado, compuesto por un valle 
con sección parabólica y un cauce con sección 
trapezoidal. En todos los casos se consideró que 
aguas abajo de la presa el canal se encuentra 
inicialmente seco. La figura 9 muestra esque- 
máticamente la geometría del sistema hidráu- 


Cuadro 4. Condiciones iniciales y de borde. 


Condición inicial seca 


Condición inicial mojada 


10, x<100 
(bx) 0, x>100 


Ux,4.) E Uy(x.10) =0 


10, x<100 
5, x>100 


(20) 


Uso) fo) 0 


200 


0 50 100 150 200 
x (m) 


Figura 8. Malla y superficie libre calculada con el modelo (izquierda) y mapa de isolíneas de profundidad (derecha), 
para canal inicialmente seco (arriba) y mojado (abajo) en t=7.2 s. 
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Figura 9. Esquema del sistema hidráulico (izquierda) e hidrograma de salida (derecha). 


lico y la forma típica del hidrograma de salida 
que causa un rompimiento de la presa. 


Hidrograma de salida 


La figura 10 muestra el caudal máximo del 
hidrograma de salida en la sección de la presa, 
calculado con el modelo versus el caudal 
correspondiente a una onda de gravedad, 
calculado como: 


Ori Esta 


donde h es la profundidad de la presa y A es la 
sección transversal de la brecha. Se muestran 
cálculos para profundidades de la presa de 45, 
60, 75, 90, 105 y 120 m, y anchos de brecha de 
400, 450, 500, 550, 600 y 650 m. 

Se observa que existe una relación lineal 
entre ambos caudales, de manera que el caudal 
máximo a la salida podría estimarse como 


Oré =0.41 - Qonda" 


Efecto de la brecha, pendiente de fondo, rugosidad y 
orientación 


La figura 11 muestra la variación del caudal 
máximo de salida correspondiente a la situa- 
ción base (brecha = 400 m; n = 0.02 sm”; 


orientación = 0%) con la brecha, rugosidad 
del cauce y orientación del muro de la presa 
respecto al eje longitudinal del valle. 

Se observa que la rugosidad y orientación 
tienen un efecto menor en el caudal máximo del 
hidrograma de salida. Sin embargo, el caudal 
máximo crece en forma rápida con el ancho de 
la brecha hasta que se produce el rompimiento 
total. 


Atenuación del hidrograma de salida 


La atenuación del hidrograma de salida hacia 
aguas abajo se estudió mediante la simulación 
de 28 escenarios, en que se varió la profundidad 
del embalse, orientación de la presa, ancho de 
la brecha y pendiente longitudinal respecto al 
eje longitudinal del valle, a fin de analizar la 
sensibilidad de la respuesta hidráulica a las 
propiedades geométricas del sistema. En cada 
caso se cuantificó el tiempo base y el caudal 
máximo del hidrograma, en un punto ubicado 
2 km aguas abajo de la presa. El cuadro 5 
muestra los escenarios analizados y la figura 12 
muestra el caudal máximo y tiempo de crecida 
para dichos escenarios. 

Se observa, al igual que en el hidrograma 
de salida en la sección de la presa, que la 
variable más sensible fue la profundidad de 
la presa, que al variar de 45 a 120 m provocó 
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Figura 10. Caudal máximo del hidrograma y caudal de una onda de gravedad. 
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Figura 11. Efecto de la brecha (izquierda), rugosidad (centro), y orientación de la presa (derecha) sobre el caudal máximo. 


cambios en el caudal máximo de 20 000 a 90 000 
m?*/s, seguida de la brecha que en el caso de 
rompimientos parciales puede hacer variar el 
caudal máximo en un factor tres. La rugosidad 
del cauce puede disminuir 
considerable el caudal máximo y aumentar 


de manera 
el tiempo de crecida. La pendiente de fondo 
y Orientación tienen un efecto menor sobre el 
hidrograma y su atenuación hacia aguas abajo. 


Conclusiones 


Se aplicó el método de los volúmenes finitos 
para dar solución a las ecuaciones de Saint 


Venant 2D con términos fuente, a fin de calcular 
el flujo inducido por el rompimiento de una 
presa. 

El modelo se verificó mediante el cálculo 
de siete casos con solución conocida, que 
permitieron determinar la capacidad para la 
simulación de escurrimientos con disconti- 
nuidades y simulación del proceso secado / 
mojado. 

Se analizó el flujo inducido por el rom- 
pimiento instantáneo de una presa en un 
sistema hidráulico idealizado, compuesto por 
un valle con sección parabólica y un cauce con 
sección trapezoidal. Se simularon 28 escenarios, 
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Cuadro 5. Escenarios analizados. 


Escenario Profundidad Pendiente Brecha Orientación Rugosidad 
Escenario 01* 45 0.005 100 0 0.025 
Escenario 02 45 0.005 100 15 0.025 
Escenario 03 45 0.005 100 30 0.025 
Escenario 04 45 0.005 100 45 0.025 
Escenario 05 45 0.005 100 60 0.025 
Escenario 06 45 0.000 100 0 0.025 
Escenario 07 45 0.0025 100 0 0.025 
Escenario 08 60 0.005 100 0 0.025 
Escenario 09 75 0.005 100 0 0.025 
Escenario 10 90 0.005 100 0 0.025 
Escenario 11 105 0.005 100 0 0.025 
Escenario 12 120 0.005 100 0 0.025 
Escenario 13 45 0.005 50 0 0.025 
Escenario 14 45 0.005 100 0 0.025 
Escenario 15 45 0.005 200 0 0.025 
Escenario 16 45 0.005 300 0 0.025 
Escenario 17 45 0.005 400 0 0.025 
Escenario 18 45 0.005 100 0 0.000 
Escenario 19 45 0.005 100 0 0.010 
Escenario 20 45 0.005 100 0 0.020 
Escenario 21 45 0.005 100 0 0.030 
Escenario 22 45 0.005 100 0 0.040 
Escenario 23 45 0.005 100 0 0.050 
Escenario 24 45 0.005 100 0 0.060 
Escenario 25 45 0.001 100 0 0.025 
Escenario 26 45 0.002 100 0 0.025 
Escenario 27 45 0.003 100 0 0.025 
Escenario 28 45 0.004 100 0 0.025 


*Caso base. 


donde se varió la altura y orientación del muro 
de la presa, el ancho de la brecha, la pendiente 
longitudinal y la rugosidad del cauce aguas 
abajo. Los resultados muestran que la altura 
de la presa es la variable más influyente en el 
caudal máximo del hidrograma de salida y que 
el ancho de la brecha tiene un efecto importante 
en rompimientos parciales. La rugosidad 
del lecho atenúa de manera considerable el 
caudal máximo hacia aguas abajo y aumenta el 
tiempo de crecida. La orientación del muro y la 


pendiente longitudinal tienen un efecto menor 
sobre las crecidas provocadas por rompimiento. 

Los resultados obtenidos pueden aplicarse 
en la planificación de presas y la gestión del 
riesgo de inundación ante rompimientos. 
El modelo desarrollado calculó de forma 
adecuada el complejo comportamiento que 
exhibe el flujo inducido por el rompimiento 
instantáneo de una presa, constituyendo una 
herramienta de simulación con capacidad 
predictiva. Como dato de entrada para el 
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Figura 12. Caudal máximo (izquierda) y tiempo de crecida (derecha) para los escenarios analizados. 


modelo desarrollado, casos particulares de graduales, se requiere conocer la forma de la 
estudio requieren la topografía del cauce y brecha en función del tiempo. 

planicies de inundación, geometría del muro El desarrollo futuro se orientará a extender 
y rugosidad. En caso de aberturas de brecha el modelo para simular el escurrimiento en 
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forma acoplada con la evolución del lecho, 
de manera que se tomen en cuenta efectos del 
transporte de sedimentos. 
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